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Abstract 

Quantum theory imposes a strict limit on the strength of non-local correlations. It only allows 
\ for a violation of the CHSH inequality up to the value 2y/2, known as Tsirelson's bound. In this 

paper, we consider generalized CHSH inequalities based on many measurement settings with two 
' possible measurement outcomes each. We demonstrate how to prove Tsirelson bounds for any 

such generalized CHSH inequality using semidefinite programming. As an example, we show 
that for any shared entangled state and observables X\ , . . . , X n and Y\ , . . . , Y n with eigenvalues 
±1 we have | (X 1 Y 1 ) + (X 2 Y 1 ) + (X 2 Y 2 ) + (X 3 Y 2 ) + ... + (X n Y n ) - (X x Y n ) \ < 2n cos (n/{2n)) . It 
is well known that there exist observables such that equality can be achieved. However, we show 
t-H . that these are indeed optimal. Our approach can easily be generalized to other inequalities for 

' such observables. 

^ : 

Non-local correlations arise as the result of measurements performed on a quantum system 
shared between two spatially separated parties. Imagine two parties, Alice and Bob, who are given 
access to a shared quantum state |^), but cannot communicate. In the simplest case, each of them 
is able to perform one of two possible measurements. Every measurement has two possible outcomes 
labeled ±1. Alice and Bob now measure 1^) using an independently chosen measurement setting 
and record their outcomes. In order to obtain an accurate estimate for the correlation between 
their measurement settings and the measurement outcomes, they perform this experiment many 
^ . times using an identically prepared state \^) in each round. Both classical and quantum theories 

impose limits on the strength of such non-local correlations. In particular, both do not violate the 
non-signaling condition of special relativity. That is, the local choice of measurement setting does 
not allow Alice and Bob to transmit information. Limits on the strength of correlations which are 
possible in the framework of any classical theory, i.e. a framework based on local hidden variables, 
are known as Bell inequalities 1 . The best known Bell inequality is the Clauser, Home, Shimony 
and Holt (CHSH) inequality ^ 

KXxiq) + {X X Y 2 ) + (X 2 Y 1 ) - (X 2 Y 2 )\ < 2, 

where {Xi,X 2 } and {Yi,Y"2} are the observables representing the measurement settings of Alice 
and Bob respectively. {XiYj) = (^>\Xi ® ij|vl/) denotes the mean value of Xi and Yj. Quantum 
mechanics allows for a violation of the CHSH inequality, but curiously still limits the strength of 
nonlocal correlations. Tsirelson's bound ^7] says that for quantum mechanics 

HXxYl) + (X{Y 2 ) + {X 2 Y{) - (X 2 Y 2 )\ < 2V2. 
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Peres demonstrated how to derive Bell inequalities ^2] even for more than two settings. As 
Froissart and Tsirelson JS] have shown, these inequalities correspond to the faces of a polytope. 
Computing the boundary of the space of correlations that can be attained using a classical theory 
therefore corresponds to determining the faces of this polytope. However, determining bounds 
on the correlations that quantum theory allows remains an even more difficult problem 4 . All 
Tsirelson's bounds are known for CHSH-type inequalities (also known as correlation inequalities) 
with two measurement settings and two outcomes for both Alice and Bob dTJ] . Filipp and Svozil [7] 
have considered the case of three measurement settings analytically and conducted numerical studies 
for a larger number of settings. Finally, Buhrman and Massar have shown a bound for a generalized 
CHSH inequality using three measurement settings with three outcomes each 

In this paper, we investigate the case where Alice and Bob can choose from n measurement 
settings with two outcomes each. We use a completely different approach based on semidefinite 
programming in combination with Tsirelson's seminal results |1T | llf>| fTTTj . This method is similar to 
methods used in computer science for the two-way partitioning problem j2] and the approximation 
algorithm for MAXCUT by Goemans and Williamson 9... Cleve et al. |Hj have also remarked 
that Tsirelson's constructions leads to an approach by semidefinite programming in the context 
of multiple interactive proof systems with entanglement. Semidefinite programming allows for 
an efficient way to approximate Tsirelson's bounds for any CHSH-type inequalities numerically. 
However, it can also be used to prove Tsirelson type bounds analytically. As an illustration, we first 
give an alternative proof of Tsirelson's original bound using semidefinite programming. We then 
prove a new Tsirelson's bound for the following generalized CHSH inequality 13]. Classically, it 
can be shown that 

n n—1 

I ^(XM) + ^(Xi+iYi) - (X 1 Y n )\ < 2n - 2. 

i=\ i=l 

Here, we show that for quantum mechanics 

n n—1 

| ^{Xtf) + 2<*i+i*i> - (XiY n )\ < 2ncos (^-) , 

i=l i=l 

where {X\, . . . , X n } and {Y"i, . . . ,Y n } are observables with eigenvalues ±1 employed by Alice and 
Bob respectively, corresponding to their n possible measurement settings. It is well known that 
this bound can be achieved |3j for a specific set of measurement settings if Alice and Bob share 
a singlet state. Here, we show that this bound is indeed optimal for any state |^) and choice 
of measurement settings. This method generalizes to other CHSH inequalities, for example, the 
inequality considered by Gisin 8 . As outlined below, Tsirelson's results also imply that any bound 
proved using this method can indeed be achieved using quantum mechanics. As Braunstein and 
Caves J3J have shown, it is interesting to consider inequalities based on many measurement settings, 
in particular, the chained CHSH inequality above. The gap between the classical and the quantum 
bound for this inequality is larger than for the original CHSH inequality with only two measurement 
settings. They show that even for real experiments that inevitably include noise, this inequality 
leads to a stronger violation of local realism, and may thus lead to a better test. 
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1 Preliminaries 



Throughout this paper, we write u = (ui, . . . , u n ) for an n-element vector, u ■ v denotes the inner 
product between vectors u and v. Furthermore, diag(A) denotes the matrix with the components 
of the vector A on its diagonal. We write A = [aij] to indicate that the entry in the i-th row and 
j-th. column of A is a^. We also use the shorthand [n] = {1, . . . , n}. A^ is the conjugate transpose 
of matrix A. A positive semidefinite n x n matrix A is a nonsingular Hermitian matrix such that 
x* Ax > for all x G C n |10|. We use A y to indicate that A is positive semidefinite. It will be 
important that a Hermitian matrix A is positive semidefinite if and only if all of its eigenvalues are 
nonnegative [TUJ Theorem 7.2.1]. 

We will need two ingredients for our proof. First, the following result by Tsirelson |17[ Theorem 
1] and El plays an essential role. 

Theorem 1 (Tsirelson) Let X\, . . . ,X n and Y\,...,Y n be observables with eigenvalues in the 
interval [—1; 1]. Then for any state |^) G A <S> B and for all s,t G [n] there exist real unit vectors 
x\, . . . , x n ,yi, . . . , y n G M 2n such that 

(^\X s ®Y t \ty) =x s -y t . 

Conversely, let x s ,yt G l w be real unit vectors. Let \^>) G A(£>B be any maximally entangled state 
where dim(^4) = dim(B) = 2^ N ^ 2 \ Then for all s,t G [n] there exist observables X s on A and Y t 
on B with eigenvalues ±1 such that 

x s ■ y t = {^\X S (8) YtYSl). 

In particular, this means that we can rewrite CHSH inequalities in terms of vectors. The second 
part of Tsirelson's result implies that any strategy based on vectors can indeed be implemented 
using quantum measurements. See ^0] for a detailed construction. 

Secondly, we will make use of semidefinite programming. This is a special case of convex 
optimization. We refer to |2j for an in-depth introduction. The goal of semidefinite programming 
is to solve he following semidefinite program (SDP) in terms of the variable X G S n 

maximize Tr(CX) 

subject to Tr(AiX) = 6j, i = 1, . . . ,p, and X y 

for given matrices C, A±, . . . , A p G S n where S n is the space of symmetric n x n matrices. X is 
called feasible, if it satisfies all constraints. An important aspect of semidefinite programming is 
duality. Intuitively, the idea behind Lagrangian duality is to extend the objective function (here 
Tr(CX)) with a weighted sum of the constraints in such a way, that we will be penalized if the 
constraints are not fulfilled. The weights then correspond to the dual variables. Optimizing over 
these weights then gives rise to the dual problem. The original problem is called the primal problem. 
An example of this approach is given in the next section. Let d' denote the optimal value of the 
dual problem, and p' the optimal value of the primal problem from above. Weak duality says that 
d! > p'. In particular, if we have d' = p' for a feasible dual and primal solution respectively, we can 
conclude that both solutions are optimal. 
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2 Tsirelson's bound 



To illustrate our approach we first give a detailed proof of Tsirelson's bound using semidefinite 
programming. This proof is more complicated than Tsirelson's original proof, however, it serves 
as a good introduction to the following section. Let Xi,X 2 and Y\,Y 2 denote the observables with 
eigenvalues ±1 used by Alice and Bob respectively. Our goal is now to show an upper bound for 

KXiYx) + (X 1 Y 2 ) + (X 2 Y 1 ) - (X 2 Y 2 )\. 

From Theorem^we know that there exist real unit vectors x s , y% £ M 4 such that for all s,t £ {0, 1} 
(X s Y t ) = x s • yt- In order to find Tsirelson's bound, we thus want to solve the following problem: 
maximize xi-yi+xi-y 2 +x 2 -yi—x 2 -y 2 , subject to || x\ \\ = \\ x 2 \\ = \\ y± \\ = \\ y 2 \\ = 1. Note that we 
can drop the absolute value since any set of vectors maximizing the above equation, simultaneously 
leads to a set of vectors minimizing it by taking —y\, —y 2 instead. We will now phrase this as a 
semidefinite program. Let G = [gij\ be the Gram matrix of the vectors {x%, x 2 , y\, y 2 } C M 4 with 
respect to the inner product: 



G = 



( 


Xi 


Xi 


X\ 


x 2 


X\ 


yi 


X\ 


V2 \ 




x 2 


X 1 


x 2 


x 2 


x 2 


yi 


x 2 


V2 




yi 


X 1 


yi 


x 2 


yi 




yi 


V2 


V 


V2 


X 1 


V2 


x 2 


V2 


2/1 


V2 


V2 ) 



G can thus be written as G = B T B where the columns of B are the vectors {x\, x 2 , yi, y 2 }. By 
Theorem 7.2.11] we can write G = B T B if and only if G is positive semidefinite. We thus impose 
the constraint that G >z 0. To make sure that we obtain unit vectors, we add the constraint that 
all diagonal entries of G must be equal to 1. Define 



W = 








1 


1 


\ 








1 


-1 




1 


1 










V i 


-1 








/ 



Note that the choice of order of the vectors in B is not unique, however, a different order only 
leads to a different W and does not change our argument. We can now rephrase our optimization 
problem as the following SDP: 

maximize ^Tr(GW) 

subject to G y and Vi,gu = 1 

We can then write for the Lagrangian 

L(G, A) = ^Tr(GW) - Tr(diag(A)(G - /)), 
where A = (Ai, A2, A3, A4). The dual function is then 

5 (A) = sup Tr (g - diag(A)) ) + Tr(diag(A)) 

J Tr(diag(A)) if \W - diag(A) < 
\ 00 otherwise 

We then obtain the following dual formulation of the SDP 
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minimize Tr(diag(A)) 
subject to ~\W + diag(A) h 

Let ft and d! denote optimal values for the primal and Lagrange dual problem respectively. From 
weak duality it follows that d! > f' . For our example, it is not difficult to see that this is 
indeed true. Let G' and A' be optimal solutions of the primal and dual problem, i.e. ft = 
^Tr(G'W) and dl = Tr(diag(A')). Recall that all entries on the diagonal of G' are 1 and thus 
Tr(diag(A')) = Tr(G'diag(A')). Then d' - p' = Tr(diag(A')) - \Ti(G'W) = Tr(G"diag(A') - 
\G'W) = Tr(G'(diag(A') — \W)) > 0, where the last inequality follows from the constraints G' >z 0, 
diag(A') - \W >: and [2J Example 2.24]. 

In order to prove Tsirelson's bound, we will now exhibit an optimal solution for both the primal 
and dual problem and then show that the value of the primal problem equals the value of the dual 
problem. The optimal solution is well known [171 1151 ITT] . Alternatively, we could easily guess the 
optimal solution based on numerical optimization by a small program for Matlab 1 and the package 
SeDuMi [2] for semidefinite programming. Consider the following solution for the primal problem 



G' = 
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1 





1 

f 


1 









1 












V2 






1 


1 


1 















V 


V2 


V2 





1 


) 



which gives rise to the primal value p' = \Tt{G'W) = 2\[2. Note that G' y since all its eigenvalues 
are nonnegative |1Q1 Theorem 7.2.1] and all its diagonal entries are 1. Thus all constraints are 
satisfied. The lower left quadrant of G' is in fact the same as the well known correlation matrix for 

2 observables |16l Equation 3.16]. Next, consider the following solution for the dual problem 

A' = -^(l, 1,1,1). 

The dual value is then d! = Tr(diag(A')) = 2^2. Because -W + diag(A') >z 0, A' satisfies the 
constraint. Since p' = d' , G' and A' are in fact optimal solutions for the primal and dual respectively. 
We can thus conclude that 

\{X 1 Y 1 ) + (X 1 Y 2 ) + (X 2 Y 1 ) - (X 2 Y 2 )\ < 2V2, 

which is Tsirelson's bound ^7]. By Theorem^ this bound is achievable. 

3 Tsirelson's bounds for more than 2 observables 

We now show how to obtain bounds for inequalities based on more than 2 observables for both Alice 
and Bob. In particular, we will prove a bound for the chained CHSH inequality for the quantum 
case. It is well known that it is possible to choose observables X%, . . . , X n and Y±, . . . ,Y n such 
that 

n n—1 

I ^(X^> + J2{X i+1 Yi) - (X 1 Y n )\ = 2ncos (^) . 

i=l i=l 

1 See http://www.cwi.nl/~wehner/tsirel/ for the Matlab example code. 
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We now show that this is optimal. Our proof is similar to the last section. However, it is more 
difficult to show feasibility for all n. 

Theorem 2 Let p G A (8) B be an arbitrary state, where A and B denote the Hilbert spaces of 
Alice and Bob. Let Xi, . . . ,X n and Y\, . . . ,Y n be observables with eigenvalues ±1 on A and B 
respectively. Then 



n 

£ 

i=i 



n—l 



(XiYi) + J2{X i+1 YJ - {X 1 Y n )\ < 2ncos 



1=1 



\2nJ 



Proof. By Theorem^ our goal is to find the maximum value for x\ ■ y\ + x% ■ y\ + X2 ■ 2/2 + ^3 • 2/2 + 



• • • + x n ■ y n - xi ■ y n , for real unit vectors xi, . . . , x n , y 1 , . . . , y n £ 



d2/i 



As above we can drop the 



absolute value. Let G = [gij] be the Gram matrix of the vectors {xi, . . . , x n , yi, . . . , y n } C R 2n . As 
before, we can thus write G = B T B, where the columns of B are the vectors {x\, . . . ,x n , yi, . . . , y n }, 
if and only if G >z 0. To ensure we obtain unit vectors, we again demand that all diagonal entries 
of G equal 1. Define n x n matrix A and 2n x 2n matrix W by 



/ 1 





A 



\ 





V-i . 

We can now phrase our maximization problem as the following SDP: 



'•• 
1 1 
1/ 



W 



At 
A 



maximize ^Tr(GW) 

subject to G >z and Vi,gu = 1 

Analog to the previous section, the dual SDP is then: 

minimize Tr(diag(A)) 
subject to -\W + diag(A) h 

Let p' and d' denote optimal values for the primal and dual problem respectively. As before, d' > p'. 



Primal We will now show that the vectors suggested in are optimal. For k € [n], choose unit 
vectors Xfci 2/fc £ K 2n to be of the form 

x k = (cos(0 fc ),sin(0 fc ),O, . . . ,0), 
y k = (cos(^jfc),sin(^fc),0, ...,0), 

where = ^(2^ — 2) and ^fc = <£r(2A; — 1). The angle between and 2/fc is given by ipk — <j>k = ^ 
and thus Xk-yu = cos (^-) . The angle between Xfc+i and 2/fc is — = and thus x^+i • yk = 



cos 



^J. Finally, the angle between —x\ and y n is 7r — -0 n = 
value of our primal problem is thus given by 

n n—l 

p' = y^Xk ■ yk + ^fc+i -yt- xi-y n 



2n 



and so — xi • y n = cos 



2n/ 



The 



fe=l 



fc=i 



2ncos ( — 1 
V2n/ 
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Let G' be the Gram matrix constructed from all vectors Xk,yk as described earlier. Note that 
our constraints are satisfied: \/i : gu = 1 and G' y 0, because G' is symmetric and of the form 

a = b t b. 

Dual Now consider the n-dimensional vector 

* = «»(£) (i,...,i). 

In order to show that this is a feasible solution to the dual problem, we have to prove that -\W + 
diag(A') y and thus the constraint is satisfied. To this end, we first show that 

Claim 1 The eigenvalues of A are given by 7 S = 1 + e ,7I "( 2s + 1 )/ n with s = 0, . . . , n — 1. 

Proof. Note that if the lower left corner of A were 1, A would be a circulant matrix ^3], i.e. each 
row of A is constructed by taking the previous row and shifting it one place to the right. We can 
use ideas from circulant matrices to guess eigenvalues 7 S with eigenvectors 



ii — In 11 ' 1 n n ~ 2 o n ~ 3 n 1 o°) 



where p s = e i7r ( 2s + 1 )/ n and s = 0, . . . , n — 1. By definition, u = (u\,U2, . . . ,u n ) is an eigenvector 
of A with eigenvalue 7 if and only if Au = ju. Here, Au = "fu if and only if 

(i) Vj G {1, . . . , n - 1} : Uj + u j+1 = jUj, 
(ii) -ui +u n = 7« n . 

Since for any j £ {1, . . . , n — 1} 

Uj +u j+1 = p n -i + p n -i- l = 



e 



-i(n-j)w(2s+l)/n^ _|_ e m(2s+l)/n \ 



(i) is satisfied. Furthermore (ii) is satisfied, since 



-til + u n = -p™ 1 + p° s 



_ e -iw(2s+l) e iw(2s+l) /n _j_ 
1 _|_ e i7r(2s+l)/n _ 

lsP° s = lsU n - 



□ 



Claim 2 The largest eigenvalue of W is given by 7 = 2 cos (^) . 

Proof. By |10[ Theorem 7.3.7], the eigenvalues of W are given by the singular values of A and 
their negatives. It follows from Claim ^ that the singular values of A are 



, /vr(2s + l)\ 

^ = \/7^7j=2 + 2cosf^— J -j 
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Considering the shape of the cosine function, it is easy to see that the largest singular value of 
A is given by 2 + 2cos(-7r/n) = 4cos 2 (7r/(2n)), the largest eigenvalue of W is \J2 + 2 cos(7r/n) = 
2cos(vr/(2n)). □ 



Since — ^ 



7}W and diag(A') are both Hermitian, Weyl's theorem jlUl Theorem 4.3.1] implies that 




where "imin{M) is the smallest eigenvalue of a matrix M. It then follows from the fact that diag(A') 
is diagonal and Claim [2] that 



Because p' = d! , G' and A' are optimal solutions for the primal and dual respectively, which 



Note that for the primal problem we are effectively dealing with 2-dimensional vectors, Xk,Uk- It 
therefore follows from Tsirelson's construction ^S] that given an EPR pair we can find observables 
such that the bound is tight. In fact, these vectors just determine the measurement directions as 
given in 

4 Discussion 

Our approach can be generalized to other CHSH-type inequalities. For another inequality, we 
merely use a different matrix A in W. For example, for Gisin's CHSH inequality [S], A is the matrix 
with l's in the upper left half and on the diagonal, and -l's in the lower right part. Otherwise our 
approach stays exactly the same, and thus we do not consider this case here. Numerical results 
provided by our Matlab example code suggest that Gisin's observables are optimal. Given the 
framework of semidefinite programming, the only difficulty in proving bounds for other inequalities 
is to determine the eigenvalues of the corresponding A, a simple matrix. Since our approach is 
based on Tsirelson's vector construction, it is limited to CHSH-type inequalities for two parties. 
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□ 



8 



Massar for pointers, to Ronald de Wolf for the pointer to |13j and proofreading, and to Manuel 
Ballester for alerting me to correct the journal reference of f3J. Supported by EU project RESQ 
IST-2001-37559 and NWO vici project 2004-2009. 



References 

[1] J. S. Bell. On the Einstein-Podolsky-Rosen paradox. Physics, 1:195-200, 1965. 

[2] S. Boyd and L. Vandenberghe. Convex Optimization. Cambridge University Press, 2004. 

[3] S.L. Braunstein and CM. Caves. Wringing out better Bell inequalities. Annals of Physics, 
202:22-56, 1990. 



[4] H. Buhrman and S. Massar. Causality and Cirel'son bounds, quant-ph/0409066, 2004. 

[5] J. Clauser, M. Home, A. Shimony, and R. Holt. Proposed experiment to test local hidden- 
variable theories. Physical Review Letters, 23:880-884, 1969. 

[6] R. Cleve, P. H0yer, B. Toner, and J. Watrous. Consequences and limits of nonlocal strategies. 
In Proceedings of 19th IEEE Conference on Computational Complexity, pages 236-249, 2004. 
|quant-ph/0404076[ 

[7] S. Filipp and K. Svozil. Tracing the bounds on Bell- type inequalities. In Proceedings of 
Foundations of Probability and Physics-3, pages 87-94, 2004. 

[8] N. Gisin. Bell inequality for arbitrary many settings of the analyzers. Physics Letters A, 
260:1-3, 1999. 

[9] M.X. Goemans and D.P. Williamson. Improved approximation algorithms for maximum 
cut and satisfiability problems using semidefinite programming. J. Assoc. Comput. Mach., 
42:1115-1145, 1995. 

[10] R. A. Horn and C. R. Johnson. Matrix Analysis. Cambridge University Press, 1985. 
[11] A. Peres. Quantum Theory: Concepts and Methods. Kluwer Academic Publishers, 1993. 
[12] A. Peres. All the Bell inequalities. Foundations of Physics, 29:589-614, 1999. 
[13] R.M.Gray. Toeplitz and Circulant Matrices: A review. 1971. 
[14] J. Sturm and AdvOL. SeDuMi. http://sedumi.mcmaster.ca/ 

[15] B. Tsirelson. Quantum analogues of Bell inequalities: The case of two spatially separated 
domains. Journal of Soviet Mathematics, 36:557-570, 1987. 

[16] B. Tsirelson. Some results and problems on quantum Bell-type inequalities. Hadronic Journal 
Supplement, 8(4):329-345, 1993. 

[17] B. Cirel'son (Tsirelson). Quantum generalizations of Bell's inequality. Letters in Mathematical 
Physics, 4:93-100, 1980. 



9 



